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Abstract 

We construct rotating hairy black holes in SU(2) Einstein- Yang-Mills theory. 
These stationary axially symmetric black holes are asymptotically flat. They 
possess non-trivial non-Abelian gauge fields outside their regular event horizon, 
and they carry non-Abelian electric charge. In the limit of vanishing angular 
momentum, they emerge from the neutral static spherically symmetric Einstein- 
Yang-Mills black holes, labelled by the node number of the gauge field function. 
With increasing angular momentum and mass, the non-Abelian electric charge of 
the solutions increases, but remains finite. The asymptotic expansion for these 
black hole solutions includes non-integer powers of the radial variable. 
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1 Introduction 



The unique family of stationary asymptotically flat black holes of Einstein-Maxwell 
(EM) theory comprises the rotating Kerr-Newman and Kerr black holes and the static 
Reissner-Nordstr0m and Schwarzschild black holes. EM black holes are completely 
determined by their mass, their charge and their angular momentum, i.e. EM black 
holes have "no hair" [jl], . 

The EM "no-hair" theorem does not generalize to theories with non-Abelian gauge 
fields coupled to gravity 0. The generic black hole solutions of SU(2) Einstein- Yang- 
Mills (EYM) theory possess non-trivial magnetic fields outside their regular event hori- 
zon, representing non-Abelian "hair" Besides static spherically symmetric hairy 
black holes Q], there are also static hairy black holes, which are not spherically but 
only axially symmetric |^, 0]. This shows that Israel's theorem neither generalizes 
to theories with non-Abelian gauge fields coupled to gravity. 

Obviously, also hairy stationary black hole solutions, representing the non-Abelian 
generalizations of the Kerr-Newman black hole solutions, should exist, as conjectured 
long ago [§. The construction of such hairy rotating black hole solutions, however, 
appeared very difficult. First of all, it was not clear, whether the standard Lewis- 
Papapetrou parametrization of the stationary axially symmetric metric would be 
sufficiently general and whether an ansatz for the gauge fields, satisfying the Ricci 
circularity and Frobenius conditions |10|, was available 0, |T^. Second, even 
within the standard metric parametrization the black hole solutions would only possess 
axial symmetry. Therefore the construction of such solutions would involve the solution 
of a large system of coupled non-linear partial differential equations for the metric and 
gauge field functions and thus represent a numerical challenge. 

First progress was achieved with the construction of perturbative stationary non- 
Abelian black hole solutions in SU(2) EYM theory [|12|. Based on the static hairy 
black hole solutions, these slowly rotating black hole solutions revealed an unexpected 
property. The rotating black holes carry non-Abelian electric charge, whereas their 



static counterparts are electrically neutral [12]. In the static spherically symmetric 
case electrically charged SU(2) black holes are even prohibited by the 'non-Abelian 



baldness' theorem [|14[. Indeed, the non-Abelian electric charge of the slowly rotating 
black hole solutions turned out to be proportional to their angular momentum, and 
thus vanishes in the static limit. Subsequently, perturbative calculations with more 
general boundary conditions predicted even more exotic stationary hairy black hole 
solutions p5[1 . 

Only recently non-pertubative rotating hairy black hole solutions were obtained 
in SU(2) EYM theory []TB[, confirming the perturbative calculations [0. They are 
obtained within the standard Lewis-Papapetrou parametrization of the metric, and the 
ansatz for the gauge fields is consistent with the circularity and Frobenius conditions. 
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Representing the first set of non-perturbative rotating liairy black liole solutions, they 
possess three global charges, a mass, an angular momentum, and a small non-Abelian 
electric charge. They do not carry non-Abelian magnetic charge, although they possess 
non-trivial magnetic gauge fields outside their regular event horizon. 

Here we present a detailed account of these rotating hairy black holes, announced 
in [|16l. We analyze their properties. In particular we discuss their global charges 



and their horizon charges |]T3- We further introduce local charges, to illustrate the 



contributions of the gauge fields outside the event horizon. The black hole solutions 
depend on two continuous parameters, the horizon size, and the angular velocity of the 
horizon. Like their static spherically symmetric counterparts, they further depend on 
the integer node number k of the gauge field functions. Thus for a given horizon size, 
when a small angular velocity of the horizon is imposed, a sequence of rotating black 
hole solutions, labelled by the node number k, emerges from the sequence of static 
black hole solutions. 

Furthermore, we present the expansions of the metric and gauge field functions at 
the horizon and at infinity for these black holes. The asymptotic expansion for these 
black holes has the surprising feature, that the magnetic gauge field functions approach 
their asymptotic values with non-integer powers of the radial coordinate. In particular, 
the non-integer powers depend on the non-Abelian electric charge. In the static limit, 
the non-Abelian charge vanishes and the well-known integer power fall-off of the static 
spherically symmetric functions is recovered. 

In section 2 we recall the SU(2) EYM action and the equations of motion. We 
present the stationary ansatz for the metric and the gauge field, and discuss its residual 
U(l) gauge invariance. The global properties of the black hole solutions and their 
horizon properties are presented in section 3. They are obtained from the expansions at 
infinity and at the horizon. The expansions also suggest the set of boundary conditions 
to be satisfied by the solutions at the horizon and at infinity. The boundary conditions 
along the axes, follow from symmetry and regularity conditions. For comparison, we 
present in section 4 the embedded Kerr-Newman solutions. In particular, we transform 
the Kerr-Newman solutions from Boyer-Lindquist to isotropic coordinates, and further 
transform them to the gauge employed for the non-Abelian solutions. Our numerical 
results are discussed in section 5. In section 6 we present our conclusions. Appendix A 
demonstrates the circularity condition. Appendix B and C give details of the expansion 
at infinity and at the origin, respectively. 



2 SU(2) EYM action and stationary Ansatz 

We here briefly recall the SU(2) EYM action and the general set of EYM equations to 
be satisfied by the stationary hairy black hole solutions. We then discuss the ansatz 
for the metric and the gauge field functions The metric chosen is the stationary 
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axially symmetric Lewis-Papapetrou metric |^ in isotropic coordinates. The ansatz 
for the gauge field represents a generahzation of the previously employed static axially 
symmetric parameterization |T^, |5| to the stationary case, satisfying the Ricci circularity 
and Frobenius conditions 0. 

2.1 SU(2) EYM equations 

We consider the SU(2) Einstein- Yang-Mills action 

S = I (t^ - lMF,.Fn) V^d^x (1) 



Vl67rG 2 

with curvature scalar R, Newton's constant G, and SU(2) field strength tensor 

Ff,u = d^A^ - d^A^ + %e [A^, A^\ , (2) 
where A^J_ denotes the gauge fields 

K = Ir^A", (3) 
and e the Yang-Mills coupling constant. The gauge fields transform as 

A'^ = UA,U^+'-{d^U)U^ (4) 

under SU(2) gauge transformations U. 

Variation of the action ([^) with respect to the metric g'^'^ leads to the Einstein 
equations 

G^u = Rfiu — -^9tJ.vR = SnGTfj^i, (5) 

with stress-energy tensor 

dLM 



dg 



= 2Tr{F,^F,pg^^ - -^g^.F^pF'^^) , (6) 

variation with respect to the gauge field A^ leads to the matter field equations, 

1 



-g 



D.i^F^ = . (7) 
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2.2 Stationary ansatz for the metric 

To construct rotating axially symmetric EYM black hole solutions, we employ isotropic 
coordinates for the metric. In terms of spherical coordinates r, 9 and if the Lewis- 
Papapetrou metric is parameterized as ||T^ 

i,^ = -fit' + + ^dS^ + (d^ + ^dt) " , (8) 

where the four metric functions /, m, / and u depend only on the coordinates r and 9. 
The 2;-axis represents the symmetry axis. 

The metric has Killing vector fields ^ = dt and rj = d^p. The regularity condition 
along the 2;-axis 0, 

^if^ — 1, ^ = ^% , (9) 

requires 

m\e=o = l\0=o . (10) 
The ansatz for the metric, Eq. (||), satisfies the Ricci circularity conditions ||, |10|, 

i^R^iaipT]^] = = r]^Rf,io.^pr]^] , (11) 
where R^u is the Ricci tensor, and the Frobenius conditions 

^[t,VuVX];r = = ri[f,^u^x];r , (12) 

implying the corresponding conditions for the stress-energy tensor for solutions of the 
EYM equations (see section 2.3.2 and Appendix A). 

The event horizon of stationary black hole solutions resides at a surface of constant 
radial coordinate, r = rn, and is characterized by the condition /(rn) = |jl6[. The 
Killing vector field 

X = dt - uju/rud^ (13) 

is orthogonal to and null on the horizon ||19|| . 

The ergosphere, defined as the region in which ^^^^ is positive, is bounded by the 
event horizon and by the surface where 

- / + sin^ 9ju^ = . (14) 
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2.3 Stationary ansatz for the gauge field 



For the gauge field we choose the ansatz 



Af.dx^' = "^dt + A^{d^ + -dt) + (—dr + (1 - H2)de) ^ 



with 



and 



^2e 2e 



sin^ 



^2e ^ '2e 



(15) 

(16) 
(17) 



Here the symbols Xp, tq and r<^ denote the dot products of the Cartesian vector of Pauli 
matrices, f = {Tx,Ty,Tz), with the spherical spatial unit vectors, 

Cr = (sin 9 cos ip, sin 9 sin ip, cos ^) , 

60 = (cos 6* cos ip, cos 6* sin y9, — sin 9) , 

= (— sin cosy?, 0) . (18) 

The two electric gauge field functions Bi and the four magnetic gauge field functions 
Hi depend only on the coordinates r and 9. 

In the static limit, the ansatz reduces to the static spherically symmetric SU(2) 
EYM ansatz, where 



l = m , uj = , H2 = H4 , Hi = H^ = Bi = 82 = ^ , 
and all non-trivial functions depend only on the radial coordinate r. 

2.3.1 Residual gauge invar iance 



(19) 



The ansatz (|T5])-(|T^) is axially symmetric in the sense, that a rotation around the 
symmetry axis can be compensated by a gauge rotation. The ansatz is form-invariant 
under Abelian gauge transformations where 0, |, 0, |l^ 



^ = exp ( -T^T{r,9) 



(20) 



The functions Hi and H2 transform inhomogeneously under such gauge transforma- 
tions, 



Hi 
H2 



Hi - rdrT , 
H2 + deV , 



(21) 
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like a 2- dimensional gauge field. The functions H3 and H4 combine to form a scalar 
doublet, {H3 + cot9, —H4), 

H3 + cote cos r(i/3 + cot^) + sin T{-Hi) , 

-Hi cos r(-if4) -sin r(iJ3 + cote) . (22) 

Similarly, the functions Bi and B2 transform as 

-Bi + cos^— cosr(_Bi + cos^— ) + sinr(_B2 — sin6'— ) , 

ry ry ry 

Bi-sinO- ^ cosTf^a-sin^-) -sinr(5i + cose-) . (23) 
r r r 

As previously |]18|, H, 0], we choose the gauge condition 

rdrHi - deH2 = (24) 
with respect to this residual gauge degree of freedom. 

2.3.2 Stress-energy tensor 

Let us now address the stress-energy tensor T^,^ and the consistency of the ansatz for 
the gauge fields, Eqs. (p!5D-(p!7D, with the Ricci circularity conditions. The stress-energy 
is circular, when 0] 

eT.i^^m = = v'^T.^^^m ■ (25) 

To verify these circularity conditions for the stress-energy tensor T^^,, Eq. (|]), we 
expand the field strength tensor in the form 

i^.. = E^/i"^§ ' ^ = r,e,^, (26) 
where its non-vanishing components Fj^"^ are given by 
F^P = -l[H,,e + rH2,r] , 

= - sine [H,^e + H.cote + - 1] , 

F^'J = sme[H,,e + cote{H^-H2)-H2Hs] , 
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tr 



■ te 

tip 



UJ 



rBi^r + H1B2 sin e{Hi{l - H^) - H3 + rH^^r) - sin OH-i 



UJ 



rB2 r - HiBi + - sin eiHiH^ + (1 - HA +rHir)-uJr sin ^(1 - HA 
r 

B^^e - H2B2 + ^ sin ^(1/2(1 - HA) - cot9Hs - H^fi) - ^ sin^i/g 
B2,e + H2B1 - ^ sin ^(/fa^s + cot0(l - HA) - H^^e) - ^ sin Oil - H^) 



— sin 9 



UJ 



BiHi + Esfi/s + cot^) - - sin 0(cot^(l - if4) + ^^"3) 

r 



(27) 



Now it is easily seen, that Tt^ = Tte = T^r = T^e = 0. Thus the circularity conditions 
are satisfied (see Appendix A). 

Of particular interest are also the energy density of the matter fields, e = —Tq, and 
the angular momentum density j = T°. They are given by 



rpO 



and 



1 



2e^r'^m | m 



2 / _f^\\2 



UJ 



/ sin' e 



UJ 

+ - 



te 



(29) 



respectively. As seen from Eqs. (^) and (^), regularity of the energy density on the 
2;-axis requires 

H2\e=o = -f^4|6»=o • (30) 



3 Black Hole properties 

To obtain stationary axially symmetric black hole solutions which are asymptotically 
flat, and possess a regular event horizon, as well as a flnite mass, angular momentum 
and electric charge, we need to impose the appropriate set of boundary conditions. 
Looking for black hole solutions with parity reflection symmetry, we need to consider 
the solutions only for < 6 < n /2. Thus appropriate boundary conditions must be 
imposed at inflnity and at the horizon, along the p-axis and along the z-axis (i. e. for 
e = 7r/2 and 6 = 0). 

In the following we present these boundary conditions at inflnity, at the horizon, 
and along the axes. The choice of boundary conditions is consistent with the equations 
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of motion, as seen from the expansions of the metric and gauge field functions at 
infinity and at the horizon. These expansions also allow the extraction of the physical 
properties of these solutions, such as their global charges or their horizon charges. We 
also introduce local charges of the black hole solutions. 

3.1 Behaviour at infinity 

3.1.1 Boundary conditions at infinity 

For notational simplicity we now introduce the dimensionless coordinate x, 

and the dimensionless electric gauge field functions i?i and i?2, 

Bi = Bi , B2 = B2 . (32) 

e e 

To obtain asymptotically flat solutions, we impose on the metric functions at infinity 
{x = 00) the boundary conditions 

f\x=oo ^\x=oo ^|a;=oo 1 ; ^|a;=oo . (33) 

By requiring the four magnetic gauge field functions Hi to satisfy 

Hl\x=oo = -f^sU-oo = , H2\x=oo = -f^4U'=oo = =tl ; (34) 

and the two electric gauge field functions Bi to satisfy 

Bl\x=oo = B2\x=oo = , (35) 

the black hole solutions are magnetically neutral, but they may carry a non-Abelian 
electric charge. 

The node number k of the gauge field functions is defined by the number of nodes 
of the functions H2 and [|18], |^. For each node number there are two degenerate 



solutions, specified by the value of H2 and H4 at infinity, H2{oo) = -^4(00) = ±1. 
These solutions are related by the large gauge transformation U = Tr, transforming the 
gauge field functions according to 

Hi —Hi , H2 —H2 , H^ +H3 , H4 —H4 
Bi +Bi , B2 -B2 + sinO- 

X 

Because of this gauge symmetry we can choose the gauge field functions to satisfy 
-^^2(00) = H^^oo) = —1. Solutions with an odd number of nodes then have positive 
values of H2 and at the horizon, whereas solutions with an even number of nodes 
have negative values of H2 and if4 at the horizon. 
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^2 n/r2 



m — — - H sm i) -\- o 



x^ 



2aM ^ ( 1 
oj = ^ + 



x^ 



x^ 



x-" 



Hi 



2C5 , 8C4 -1(13-1) _ 1C2C-i[a + 3) ^_i(^_^^) 



3.1.2 Expansion at infinity 

The asymptotic expressions for the metric and gauge field functions are 

+ (36) 

(37) 
(38) 
(39) 

(40) 
(41) 

(42) 
(43) 
(44) 

X Vx"/ ^"^^^ 

where a, M, Q, and Ci, ... ,6*5 are dimensionless constants, and a and /3 determine 
the non- integer fall-off of the magnetic gauge field functions with 



x^ P-1 (a + 5)g2 

//2 = -l + C3X-^("-^) + o(x-5("-i)) , 

i/s = + C3X-^(°-^^) sin ^ cos ^ + o Q 
V X 



sin ^ cos 9 + 



x^ 







) sin^ ^ 









a = p-AQ\ (5=ph- AQ^ . (46) 
Further details of the asymptotic expansion are presented in Appendix B. 
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3.1.3 Global charges 

Let us now obtain the global charges of the EYM solutions from the asymptotic be- 
haviour of the metric and gauge field functions. 

The mass and the angular momentum are obtained from the metric components Qu 
and gt,p, respectively ||2^. The asymptotic expansions for the metric functions / and u, 
Eqs. (^Uj) and (P^), yield for the dimensionless mass M and the dimensionless angular 
momentum J = aM the expressions 

1 1 
M = - lim x^dx f , J = - lim x^uj . (47) 

The asymptotic behaviour of the gauge fields yields the non-Abelian electric charge 
and magnetic charge p^"^ of the black hole solutions. Let us define the gauge- 
invariant non-Abelian electric charge |F7| 




where the integral is evaluated at spatial infinity, and *F represents the dual field 
strength tensor. Insertion of the asymptotic expansion of the gauge field functions, 
Eqs. (|iO|)-(^5D, into the respective field strength tensor components, Eq. (|26|), then 
yields for the dimensionless non-Abelian electric charge 

Q= lim X f cos QB^ + sin QB^ , (49) 

i. e. the non-Abelian electric charge can be read off directly from the asymptotic be- 
haviour of the electric gauge field functions B^ and B^- 

Likewise, the gauge-invariant non-Abelian magnetic charge P^"^ is obtained from 

m 




Fl^\Qd^=^, (50) 



where the integral is evaluated at spatial infinity. Again, insertion of the asymptotic 
expansion of the gauge field functions into the respective field strength tensor compo- 
nents yields the global charge. For the dimensionless non-Abelian magnetic charge we 
thus obtain 

P = . (51) 
As imposed by the boundary conditions, the EYM solutions carry no magnetic charge. 
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3.2 Boundary conditions along the axes 

The boundary conditions along the p- and z-axis {6 = tt/2 and ^ = 0) are determined 
by the symmetries. The metric functions satisfy along the axes 

def\e=o = dem\e=o = del\e=o = deuj\e=o = , 

def\e=^ = <9em|e=| = del\e=^ = 5et^|e=f = . (52) 

For the gauge field functions symmetry considerations lead to the boundary conditions 

i^l|0=O = -f^3|e=0 = , dgH2\e=0 = dgH^\g=Q = , 

B2\g=o = , deBi\g=Q = , 

Hi\g=^ = -ff3|e=f = , deH2\e=^ = deHi\e=iL = , 

Bi\e=^=0, dgB2\e=^=0 (53) 

along the axes. In addition, regularity on the 2;-axis requires condition (0) for the 
metric functions to be satisfied, and condition ( pO| ) for the gauge field functions. 



3.3 Behaviour at the horizon 
3.3.1 Boundary conditions at the horizon 

The event horizon of stationary black hole solutions resides at a surface of constant 



radial coordinate, x = xh, and is characterized by the condition /(xh) = jlG 



Requiring the horizon to be regular, we obtain the boundary conditions at the horizon. 
At a regular horizon the metric functions must satisfy 

f\x=xii ^|x=xh ^|x=xh 5 ^|a;=a;H ) (^^) 

where is constant at the horizon. 

To obtain the boundary conditions for the gauge field functions, we start by requir- 
ing that the electro-static potential = (x'^A^), with Killing vector x, Eq. (p!3D, is 
constant at the horizon 

^'h = (X^^^)U=XH = const , (55) 
(see also Appendix C). This yields the conditions 

Bi\x=xu = -cos 6*^0 , 

B2\x=xn = sin^^^o ■ (56) 
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The equations of motion, Eqs. (0), then yield for non-Abehan solutions the conditions 

^0 = i^n/xn , (57) 
dxuj\x=xn = ^n/xH , (58) 

as well as 

{xdM\x=xu = 0, ixdxB2)U=xu = ^ . (59) 
and, taking these conditions into account, 

(deHi + xdxH2)\x=XH = 0, 
{xdxH^ - HiH4)\x=x^ = 0, 
ixdxH^ + H,{H3 + cote))\x=x^ = 0, (60) 

which derive from Fj:^, Pj:"^ , and Pj:^ , respectively. Thus the regularity conditions at 
the horizon imply 

X.^ P^i'\x=xii ) Pr6\x=xii ) PTip\x=x-ii • (6-^) 

Furthermore, the gauge condition, Eq. (^), imphes 

{xdxHi - d0H2)\x=xn = 0. 

However, for black hole solutions the gauge condition (^) still allows for non-trivial 
gauge transformations satisfying 

x^dlT + xdxT + = . (62) 

To fix the gauge, we choose the additional condition |]TB|, ^, |TB| 

ideH,)\x=xu = 0, (63) 

which implies Hi\x=xii = 0, when we take into account the boundary condition on the 
axes, i?i|e=o,7r/2 = 0, Eq. (|3|). 

For the numerical solutions we then impose on the gauge field functions the set of 
boundary conditions 

Hl\x=xu = , dxH2\x=xu = , dxH3\x=Xii = , dxH4\x=x^ = , 

xjiBi\x=xii = -COS9LJH , X}iB2\x=xji = siu 6'a;H • (64) 
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yields 



3.3.2 Expansion at the horizon 

Expanding the metric and gauge field functions at the horizon in powers of 

5 = — - 1 (65) 

Xh 

fi5,e)^5'Ml-5) + Oi5') , (66) 

m{S, 0) = 5'm2(l - 35) + 0(5^) , (67) 

l{S,e) = SH2{l-3S) + 0{S*) , (68) 

u{6,9) =u;n{l + S) + 0{6^) , (69) 

H,{5, 9)^s(l- H^, + 0{S') , (70) 

H2{S,e) = H2o + 0{6^) , (71) 
Hs{5,9) = Hso + Oi6^) , (72) 
H46,e)^Hio + 0{d^) , (73) 

5,(M) = -^^ + 0(0, (74) 

m0)^^-^ + O{5-). (75) 

The expansion coefficients /2, m2, /2, -f^ii, -f^20; -f^so; and H4Q are functions of the 
variable 9. Among these coefficients the following relations hold, 

= ^-2^, (76) 

Hu = deH2o . (77) 
Further details of the expansion at the horizon are given in Appendix C. 

3.3.3 Horizon properties 

Let us now obtain the horizon properties of the EYM solutions from the expansion at 
the horizon of the metric and gauge field functions. 

The first quantity of interest is the area of the horizon. The dimensionless area A 
is given by 

A = 27: smOy^xl . (78) 

■^0 h 
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The area A of the black hole horizon defines the area parameter xa via 



A = Attx\ . 



The entropy S of the black hole then corresponds to 



(79) 



^0) 



To obtain a measure for the deformation of the horizon we compare the dimen- 
sionless circumference of the horizon along the equator, Lg, with the dimensionless 
circumference of the horizon along the poles, Lp, 



2n 



dip 



-X sin 9 



Lp — 2 



x=xii,9=it/2 



dd 



-X 



^1) 



x=xh , V=const. 



and consider, in particular, their ratio Le/Lp 



The surface gravity of the black hole solutions is obtained from [19 



sg 



f82l 



with Killing vector x = ^ ~ (t^H/a^n)'?- Inserting the expansion in 5 = (x/xh — 1) at 
the horizon, Eqs. (|66D -(|69|), yields for the dimensionless surface gravity 



f2{0) 



(83) 



As seen from Eq. ([TUD, K^g is indeed constant on the horizon, as required by the zeroth 
law of black hole mechanics. The dimensionless temperature T of the black hole is 
proportional to the surface gravity. 



^4) 



horizon magnetic charge Pa [17|, respectively. 



Let us now consider the Yang-Mills horizon charges. By evaluating the integrals, 
Eqs. (|48|) and (pO]), at the horizon, we obtain the horizon electric charge Qa and the 



(85) 
(86) 



Q^ = i^j\Y.{*K)"ded^, 



where Qa and Pa again represent the dimensionless quantities. 
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3.4 Local charges 

To define the local mass let us first consider the general definition of the total mass 
(see e. g. ||T9[] ) 

M = 2j^ (t^, - ir^?^,) n'^i'^dV - ^ ^^e^^.^Wi'^dx^dx'' , (87) 

where S denotes an asymptotically flat hypersurface bounded by the horizon H, dV 
is the natural volume element on S, n^^ = (1, 0, 0, — u;/r)/v7 is normal to S with 



—1, and denotes the time like Killing vector field. 



Note, that T = for the action, Eq. (|T]). Straightforward calculation yields 

1 



T^un^CdV = -T^^^drdOd^ 



8nG 



Rq^/ —gdrdOdip 



where in the last step the Einstein equations have been used. 
Integration over S yields 



1 



2ti 



o n I RoV^drdOdif 
ottG je 47rG jo 



f 



r smt^ ' 



- - sin^ Ouji—-- 
dr f \drr 



dd . 



The calculation of the boundary term in Eq. 
evaluated at r = rn. 
Consequently , 



yields the same expression, but 



M 



lim — — 



/ 

1 1 /-f 



^ 2 ■ Q I J '' ■ 2 Q 

r sm — sm Uuj 



dr f 



'duo uj" 
dr r 



de 



lim , 

e G 2 Jo 

, \/4^ 1 
hm — 

e G 



Vl 2 ■ n(9f I . 2n (d^^ 

—X sm fc' sm \- 

/ \OX J \OX X . 



de 



M{x) , 



where we changed to the dimensionless coordinate x, Eq. (pi]), and defined the dimen- 
sionless local mass M{x), 



M{x) = -i -X sm^ ---sm Ou i- - - 



de 



2h f 

In a similar way we obtain the dimensionless local angular momentum J(x), 

1 



(88) 



J(x) 



7^' 
4 JO 



^9) 
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from the expression for the total angular momentum, 

1 /■ „ 1 r 1 



8nG 



L ^-"'^^'^ - 16^ /h 2 V^-'^^x^^ 



Defining the horizon mass of the black hole Ma = M(xh) and its horizon angular 
momentum Ja = J{x-e)i we obtain the relation 



Ma = 2TS - 2— Ja 



(90) 



We further define the dimensionless local non-Abelian electric and magnetic charges, 
Q{x), and P{x), 



and 



(91) 



(92) 



respectively, where the integrals are evaluated on surfaces with fixed radial coordinate 
X. Explicitly, 



Q{x) = / 
Jo 



smO—< 



— X sin^ O^^H^ + sin 9 cos 9ujil — HA — 

ox 



X sin 6 cos 0^(1 - HA + sin^ 9ujH^ - xHAcos 6Bi - sin 6B2) - 

ox 



2 / dBi dB2 
LoHi + x sin^— hcos^ 



1 2 



dx 



dx 



+ 



dui 

xsm6 cos6-^Hs 
ox 



sin^ euj{l - H^)+x sin^ ^^(1 - H4) + sin 6 cos OajH^ 

( dB OB 
xHi(sm9]3i + cos 6*^2) + x'^ cos 9— sin6'- 



dx 



dx 



2n 1/2 



, (93) 
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- sin^ + sm0 cos ^7/3(1 + ^2) + cos^ ^(^^2 - H^) + 



sin' - sin ^ cos + sin^ OH^H^ 

oO o9 



+ 



sin 9 cos 9 



MH, 



sin' 9HM + H2) + H3 - sin ^ cos 9(H2 - HA + sin 9 cos ^— ^ + 

o9 

sin' 9—^ + sin 9 cos 9H2H4 ^ d9 . 
o9 



(94) 
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4 Embedded Kerr-Newman Black Holes 



Here we briefly recall Kerr-Newman black holes, embedded in SU(2) EYM theory |22 



For better comparison with the non-Abelian rotating black holes, we first consider the 
coordinate transformation from Boyer-Lindquist coordinates to isotropic coordinates. 
Then we perform an SU(2) gauge transformation to the gauge employed for the non- 
Abelian black holes. 



4.1 Boyer-Lindquist coordinates 

In Boyer-Lindquist coordinates the metric of Kerr-Newman black holes is given by 



-— (it + a SI 

p2 V 



sm 



2 sin 6^ / , 9 , \2 P ,_9 9 ,^9 
H — [adt + p^dif) + ^dx^ + p^de^ 



where 



X + a cos U 



p- 



2 -2,2 

Pq = X + a 



(95) 



(96) 



and M denotes the black hole mass, a the angular momentum per unit mass, a = J/M, 
and C the "total charge" (see Eq. (|8|)). 

The gauge field of embedded Kerr-Newman solutions is given by 



Al^dx'' = (dt + a sin^ Od^p^ + ^"^'^ (^adt + pld^p^ , 



P" cos 9 



(97) 

y ■ / - . - . , 

where and are constant vectors in the Lie algebra, considered the Yang-Mills 
analogue of the electric and magnetic charge, respectively l^^j, ^|. They define the 
"total charge" C, 

= Q'^Q" + P'^P'^ , (98) 

and is proportional to P° |^ . 

The condition A(a;H) = yields the regular event horizon of the Kerr-Newman 
solutions, 

xh = M + Jm^ - (a2 + C2) , (99) 



extremal Kerr-Newman solutions satisfy xh = M = ^/a? + CP. 



4.2 Isotropic coordinates 

Let us now consider the Kerr-Newman solution in isotropic coordinates. For Kerr- 
Newman solutions the isotropic radial coordinate x is related to the Boyer-Lindquist 
radial coordinate x by 

, . 1 



X[X) 



£2 - 2Mx + (a2 + C2) + {x-M) 



(100) 
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and the metric functions /(x, 9), m(x, 9), l{x, 9) and uj{x, 9) (see Eq. (Q)) are given by 

f{x,e)= , , (101) 



+ cos^ 6* 



X x^ 



(102) 



and 



where 



xH(x,9) 2 2 2 ^ ^ 

cu(x,^) a(2M£-C2) 



X Pq — 0? sin ^^A 



(104) 



= ^ + M + . (105) 



In isotropic coordinates the event horizon resides at 



XH = 1/2^M2 - (a2 + C2) , (106) 

thus extremal Kerr-Newman solutions satisfy xh = 0. At the event horizon, the metric 
function / is identically zero, the metric function ratios (m/f) and (///) ^ire finite, and 
the metric function lj is constant, 

"^^^^ = 4MXH + 2M2 - = ^« ■ 
Consequently the Kerr-Newman solutions satisfy 



M2 - 4x?r - C2 



Xh 2M(M + 2xh) - (M + 2xh)2 + a- 



2 



(10^ 



4.3 Gauge transformation 



Let us now consider a Kerr-Newman solution with Lie algebra vectors Q"^ and P"^ 
pointing in the z-direction, 

(q^X / \ 7~ P COS / \ 7~ 

A^dx" = ^{dt + a sin^ 9d^) ^ + — (adt + pldip) ^ . (109) 

^ / 2 p^ ^ ^ / 



To obtain the solution in the gauge employed for the non-Abelian ansatz (|I5|)-(|I^, we 
perform a gauge transformation, where 

U{9,ip) = e-'^^^e-'^^y . (110) 
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The gauge transformed vector field reads 



A'.dx^ = + ^^^) + A^d^ + ^d9 (111) 



witfi 

A^= [^asm^e + ^p2^cos^U^-sm^^ . (112) 

\ P P ) ^ ^ 

By comparing this embedded Abelian expression for the gauge fields with the non- 
Abelian ansatz, Eqs. (|T5| ) - ([T7|) , we obtain the embedded Abelian gauge field functions 
Hi and 

C^oc P cos 
Hi = H2 = H4 = , - sin 9H3 = ■^asin'^e + ^p^ + cos^, (113) 

Bi = sm 9 — H^-l , B2 = smO — . (114) 

X X 

The asymptotic expansion for the Abelian gauge field and metric functions is given 
in Appendix B, the expansion at the horizon in Appendix C. 



5 Numerical Results 

The rotating hairy black hole solutions of SU(2) EYM theory emerge from the static 
hairy black hole solutions, when a small value of the angular velocity of the horizon is 
imposed via the boundary conditions. We therefore briefly recall the properties of the 
static hairy black hole solutions of SU(2) EYM theory Q]. 

The static hairy black hole solutions of SU(2) EYM theory carry non-trivial mag- 
netic gauge fields outside their regular event horizon, but carry no non-Abelian mag- 
netic charge. Since their electric fields vanish identically, their only global charge is 
their mass. These black holes can be characterized by their mass and by two inte- 



gers, the azimuthal winding number n of their gauge fields |]2^ and the node number 
k of their magnetic gauge field functions 0, |^, For a given winding number n, 
the solutions form sequences labelled by the node number k. With increasing k, the 
non-Abelian black hole solutions approach a limiting solution, corresponding to an 
embedded Abelian solution with magnetic charge n [|26|| . 

In contrast to the static hairy black hole solutions, their rotating generalizations 
carry electric gauge fields. Notably, the rotating hairy black hole solutions possess a 
non-Abelian electric charge |T^, [1^, whereas their static counterparts are neutral. Thus 
their global charges are their mass, their angular momentum and their non-Abelian 
electric charge. Their non-Abelian magnetic charge is identically zero. Like their static 
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counterparts, the rotating hairy black hole solutions form sequences depending on the 
winding number n and the node number k of the gauge fields. 

Here we consider the sequence of rotating hairy black hole solutions with winding 
number n = 1 and node number k. In the static limit these black hole solutions 
approach spherically symmetric SU(2) EYM black hole solutions 0]. Rotating black 
hole solutions with higher winding number n approach in the static limit SU(2) EYM 
black hole solutions, which possess only axial symmetry Such black hole solutions 
will be considered elsewhere. 

For a given node number fc, the rotating hairy black hole solutions depend on the 
isotropic horizon radius xr and on the value of the metric function uj at the horizon, 
wh, via the boundary conditions. The ratio ojr/xr represents the rotational velocity of 
the horizon. 

When the black hole solutions are constructed, their global charges are obtained 
from the asymptotic fall-off of the metric and gauge field functions of the solutions. 
Likewise, their horizon properties are obtained from the expansion at the horizon. 

To construct the rotating hairy black hole solutions, we solve the set of ten coupled 
non- linear elliptic partial differential equations numerically |2^, subject to the above 
boundary conditions. We employ compactified dimensionless coordinates, mapping 
spatial infinity to the finite value x = 1, where 

x = l-—. (115) 

X 

We furthermore introduce the functions /(x, 6), l{x, 9), rh{x, 9), and g{x, 9) [1^, |^, |^, 

f(- Q\ /(^'^) j(- n\ K^,9) -^-n^ rn{x,9) m{x,9) 
/(a;, 6) = -n , K^, d) = —— , m(x, 9) = — - — , g{x, 9) = . (116) 

x^ x^ x"^ l[x,9) 

The numerical calculations, based on the Newton-Raphson method, are performed with 



help of the program FIDISOL [|23. The equations are discretized on a non-equidistant 
grid in x and 9. Typical grids used have sizes 100 x 20, covering the integration region 
< X < 1 and < 9 < tt/2. (See ||18|, |], 0] and l2^ for further details on the numerical 
procedure.) 

5.1 Global charges 

To construct a rotating black hole solution with one node and horizon radius xh, we 
start from the static spherically symmetric SU(2) EYM black hole solution with one 
node and the same horizon radius, which represents the limiting solution in the static 
limit cuh —>■ 0. By imposing a small but finite value of cuh via the boundary conditions, 
a rotating black hole solution is obtained, which possesses non-trivial functions uj, Hi, 
H^, Bi, B2. The rotating black hole solution is only axially symmetric. It possesses a 



non-vanishing electric gauge field, giving rise to a non-Abelian electric charge |]T2|, 116 . 
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When increasing oou from zero, while keeping xu fixed, a branch of black hole 
solutions forms, the lower branch. The lower branch extends up to a maximal value 
uj^^^, where a second branch, the upper branch, bends backwards towards un = 0. 
Along both branches the mass M, the angular momentum J, and the non-Abelian 



electric charge Q continuously increase [|1^. This is seen in Figs, la-c, where the mass 



M, the angular momentum per unit mass a = J/M and the non-Abelian electric charge 
Q are shown as functions of the parameter cuh for three values of the isotropic horizon 
radius, xr = 0.1, 0.5 and 1. 

Whereas both mass M and angular momentum per unit mass a of the non-Abelian 
solutions increase strongly along the upper branch, diverging with a;^^ in the limit 
u!b_ — > 0, the non-Abelian electric charge Q remains small. It approaches apparently 
a finite limiting value, Qum ~ 0.124 (as discussed in section 5.3), independent of the 
isotropic horizon radius xh- 

For comparison, we exhibit in Figs, la-b also the mass and the angular momentum 
of embedded Abelian solutions with the same horizon radii: those of the Kerr solutions 
and those of the Kerr- Newman solutions with Q = and P = —1. (Mass and angular 
momentum of embedded Kerr-Newman solutions, possessing the same charge Q as the 
non-Abelian solutions but with P = are graphically indistinguishable from the mass 
and angular momentum of the Kerr solutions, shown. Likewise, mass and angular 
momentum of embedded Kerr-Newman solutions, possessing the same charge Q as the 
non-Abelian solutions and P = —1 are graphically indistinguishable from the mass and 
angular momentum of the Kerr-Newman solutions, shown.) 

Interestingly, both mass and angular momentum of the non-Abelian solutions, which 
carry a small electric charge Q and no magnetic charge, are close to mass and angular 
momentum of the embedded Kerr-Newman solutions with magnetic charge P = —1. 
This is seen in particular for solutions with small horizon radii xh, where large devia- 
tions of these global properties from those of the corresponding Kerr solutions arise. For 
non-Abelian black hole solutions with large horizon radii, these global charges are also 
close to those of the Kerr solutions. Here a magnetic charge of magnitude \P\ = 1 is less 
important for the global properties of the solutions. Thus also Kerr and Kerr-Newman 
solutions differ less for large horizon radii. 

As a point in case, let us inspect the maximal value of ci^h reached, u^'^^. The 
maximal value u^'^^ of the non-Abelian black hole solutions depends on the horizon 
radius xh, and increases monotonically with xh- For the solutions shown in Figs, la-c, 
u;^''^{xu = 0.1) ^ 0.0288, a;{?^^(xH = 0.5) ^ 0.0649, and u;^^'^^(xh = 1.0) ^ 0.0719. 
For Kerr-Newman solutions with Q = and P = —1 very similar values for u^^^ are 
obtained, cjf'^^(xH = 0.1) ^ 0.0278, cuf^^(xH = 0.5) ^ 0.0646, and cjf^^(xH = 1.0) ^ 
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0.0718. In contrast, for the Kerr solutions the maximal value of ujh, 



max,K 
H 



2(a/5 + 1)\ a/5 + 3 



V5 - 1 

- 0.07507 , 



is independent of the horizon radius x^. We conjecture that, in the limit xh — > oo, the 
maximal value cu^^ of the non-Abelian solutions tends to the maximal value a;^^'^ of 
the Kerr solutions. 

So far we have limited the discussion to non-Abclian solutions with one node. When 
we consider black hole solutions with higher node numbers, we observe similar features. 
The rotating hairy black hole solution with k nodes emerges from the corresponding 
static black hole solution, when a finite value of cuu is imposed via the boundary condi- 
tions. When uu increases, the lower branch of black hole solutions forms. It extends up 
to a maximal value uj^'^^, where the upper branch forms and bends backwards towards 
cjh = 0. Along both branches mass, angular momentum, and non-Abelian electric 
charge continuously increase. 

For node number A; = 3, the mass and angular momentum of the black hole solutions 
differ only little from the mass and angular momentum of the k — 1 black hole solutions. 
In general, the values are still closer to those of the Kerr-Newman solutions with Q = 
and P = —1. Also, the values of wlf^^ are very close to those of the Kerr-Newman 
solutions with Q = and P = —1. Along the lower branch, the non-Abelian electric 
charge of the k — 3 solutions is much smaller than the electric charge of the k — 1 
solutions, typically it is on the order of 10~^. Along the upper branch the non-Abehan 
electric charge increases relatively strongly. However, because of numerical inaccuracies 
encountered along the upper branch for the higher node solutions, we cannot infer from 
the numerical calculations, that the charge tends to a limiting value Qiim- 



5.2 Energy density 

Let us now discuss the rotating black hole solutions in more detail and begin with the 
energy density of the black hole solutions, e = —Tq. 

Due to the rotation, the energy density of the matter fields is angle-dependent, and, 
in particular, not constant at the horizon. The maximum of the energy density resides 
on the p-axis at the horizon, as seen in the following two representative examples. 

In Figs. 2a-d we exhibit the energy density of the matter fields of the black hole 
solution with horizon radius xa = 1 and ui-a = 0.04 on the lower branch. Fig. 2a shows 
a 3-dimensional plot of the energy density as a function of the coordinates p — xsinO 
and z = X cos 9 together with a contour plot, and Figs. 2b-d show surfaces of constant 
energy density. The surfaces of constant energy density appear ellipsoidal, being flatter 
at the poles than in the equatorial plane. For the largest values of the energy density, 
as shown in Fig. 2d, the horizon is seen in the pole region. 
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Analogously, we exhibit in Figs. 3a-d the energy density of the matter fields of the 
black hole solution with horizon radius xh = 1 and un = 0.04 on the upper branch. 
The energy density is much stronger deformed now, showing a large peak on the p-axis. 
Consequently, the surfaces of constant energy density appear torus-shaped, with the 
horizon seen in the center of the torus. Further away form the black hole horizon, 
however, the surfaces of constant energy density appear again ellipsoidal. 

We do not observe a strong dependence of the energy density on the node number 
k. On the lower branch the k = 3 densities are similar to the k = 1 densities. On the 
upper branch, the numerical determination of the k = 3 densities becomes, however, 
unreliable. 

5.3 Gauge field and metric functions 

We now turn to the metric and gauge field functions of the non-Abelian black hole 
solutions. In Figs. 4a-j we show the functions for the solutions with one node on 
the upper and lower branch for horizon radius xh = 1 and cuh = 0.04. Also shown 
are the functions of the three node solution on the lower branch for the same set of 
parameters. We first discuss the magnetic gauge field functions, shown in Figs. 4a- 
4d, then the electric gauge field functions, shown in Figs. 4e-4f, and finally the metric 
functions, shown in Figs. 4g-4j. Since the global mass and angular momentum of the 
non-Abelian solutions are rather close to those of the Kerr-Newman solutions with 
Q = and P = —1, we compare with the corresponding functions of these Kerr- 
Newman solutions. 

The magnetic gauge field function Hi, shown in Fig. 4a, is always very small. For the 
k = 1 solutions, its magnitude is two orders of magnitude smaller on the lower branch 
than on the upper branch. For the k = 3 lower branch solution it is still smaller by two 
orders of magnitude than for the lower branch k = 1 solution. Hi shows a distinct angle 
dependence, which is similar for the three black holes exhibited. The Kerr-Newman 
solutions have vanishing Hi. This is consistent with the observed smallness of Hi for 
the non-Abelian solutions. 

The magnetic gauge field function H2 is almost spherically symmetric for the k = 1 
solutions on both branches as well as for the k = 3 solution on the lower branch. It 
is exhibited in Fig. 4b. Whereas H2 decreases monotonically for the k = 1 solutions 
from a finite positive value at the horizon to its boundary value of —1 at infinity, H2 
oscillates around zero (possessing three nodes) for the k = 3 solution in an extended 
interior region, before it decreases to its boundary value —1 for x — 00. Thus H2 
precisely keeps its features of the static non-Abelian solutions, where with increasing 
node number k, it becomes close to zero in an increasing interior region. In particular, 
in the limit k —>■ 00, a limiting solution with = is reached PB| . 

The magnetic gauge field function H^, shown in Fig. 4c, is very similar for the k = 1 
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and k = 3 solutions on the lower branch, and about an order of magnitude smaller on 
the lower branch than on the upper branch. Comparison of the function of the non- 



Abelian solutions with the function H^, Eq. ( |113D , of the corresponding Kerr- Newman 
solutions with Q = and P = — 1 reveals very good agreement. Thus H^, like Hi, 
corresponds approximately to the Abelian function. 

The magnetic gauge field function H^, shown in Fig. 4d, is rather similar to the 
function H2 on the lower branch, both for the k = 1 and the k = 3 solutions. Its 
angle dependence on the lower branch is only small. On the upper branch, however, a 
distinct angle dependence appears. H4, like H2, retains the non-Abelian character of 
the solutions. 

Turning to the electric gauge field functions Bi and B2 we observe, that the angle 
dependence is largely determined by the boundary conditions at the horizon. To elim- 
inate this trivial angle dependence, we consider the new functions Bi and B2, defined 
via 

Bi = - cos dBi + sin 6132 , B2 = sin OBi + cos 9B2 . 

As seen in Fig. 4e, the function Bi then shows only a relatively small deviation from 
spherical symmetry for all three black hole solutions considered. The function B2 is 
considerably smaller in amplitude than the function Bi. Its angle dependence is similar 
for the k = 1 and k = 3 solutions on the lower branch, but different on the upper branch. 

Let us again compare with the electric gauge field functions of the corresponding 
Kerr-Newman solutions, Eq. (|114|) , with Q = and P = —1. We again find good 



agreement between these Abelian and the non-Abelian functions. Thus Bi and B2 also 
correspond in good approximation to the Abelian functions. 

Thus we conclude, that the non-Abelian gauge field functions, which are non- 
vanishing in the static limit, H2 and H4, retain their non-Abelian features. Whereas all 
gauge field functions which vanish in the static limit. Hi, H^, Bi and B2, are very close 
to the corresponding functions of the Kerr-Newman solution with Q = Q and P = — 1. 

Let us now turn to the metric functions. In order to exhibit the behaviour of the 
metric functions /, m and / at the horizon more clearly we show the functions f{x, 9), 
T{x,e), and g{x,e)_, Eq. (|lT|). 

The function / is shown in Fig. 4g. It is almost spherically symmetric on the lower 
branch and deviates from spherical symmetry only slightly on the upper branch. On the 
upper branch the magnitude of / and its slope at the horizon are considerably smaller 
than on the lower branch. There is very little dependence of the metric function / on 
the node number. This is also true for the other metric functions. 

In Fig. 4h we show the ratio g = m/l. Note, that for a spherically symmetric metric 
the function g is identical to one, thus the deviation of this function from one indicates 
the deviation from spherical symmetry. We observe that on the lower branch g is close 
to one for all x except in a region near the horizon, where it deviates slightly from one. 
On the upper branch, in contrast, the deviation from one is large and g approaches its 
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asymptotic value only for larger values of x. 

The metric function / is shown in Fig. 4i. For solutions with k = \ and A; = 3 nodes, 
as well as for solutions on the lower and upper branch, the functions / are almost 
identical. Also, their deviation from spherical symmetry is extremely small. 

We exhibit the function io in Fig. 4j. The shape of the function uj is similar for 
solutions on the lower and upper branch. However, on the upper branch the maximum 
is larger than on the lower branch. The angle dependence is small, except near the 
maximum of the function on the upper branch. 

In general, the deviation from spherical symmetry is small on the lower branch. We 
also observe little dependence of the metric functions on the node number, although 
for solutions with smaller horizon radius the node number dependence of the metric 
functions increases slightly. All metric functions are rather close to the corresponding 
functions of the Kerr- Newman solutions with (5 = and P = — 1. 

We conclude, that the functions of the non-Abelian solutions are rather close to 
those of the embedded Kerr-Newman solutions with (5 = and P = — 1, except for 
those gauge field functions, which do not vanish in the static limit. It is surprising, 
that already the one node solutions follows so closely these Abelian solutions, because 
their functions H-i and are still very different from zero, the value assumed by the 
Abelian solutions, Eq. ( |113| ). With increasing node number, though, one expects the 
solutions to get close to the Abelian solution, since and H/^ are close to zero in an 
increasing interval. 

Let us now address the limiting non-Abelian solutions, obtained in the limit cjh — > 

0. Along the lower branch the limiting solution corresponds to a static spherically 
symmetric non-Abelian black hole solution. In contrast, along the upper branch mass 
and angular momentum of the solution diverges, thus the limiting solution is singular. 
Comparison with the Abelian solutions shows, that in the limit — 0, electric and 
magnetic charge become negligible when Q ^ \ and |P| = 1. Consequently, the ratio 
Mja tends to one for fixed xh- 

We gain some understanding of the limiting solution by noting, that it can be ob- 
tained numerically, as a non-Abelian gauge field solution in an extremal background 
metric. Scaling the radial coordinate and angular momentum per unit mass by the 
mass, and taking the limit M — > oo, the metric decouples from the gauge field, 

1. e. G^J^y = 0. Taking the limiting extremal Kerr solution as background metric, one 
can solve the gauge field equations numerically. The corresponding solution is then in- 
dependent of the isotropic horizon radius xh- The limiting value of the electric charge 
obtained in this way is Qum = 0.124. 
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5.4 Local charges 



In Figs. 5a-d we present the local charges M{x), J{x), Q{x) and P{x), Eqs. 
and Eqs. (p3|)-(p^), for the black hole solutions with one node on the lower and upper 
branch, for horizon radius Xh = 1 and cjh = 0.04. We also show these charges for the 
black hole solution with three nodes on the lower branch for the same set of parameters. 

In Fig. 5a the local mass M{x) is shown. We observe, that the mass at infinity, 
representing the global mass, is only slightly larger than the mass at the horizon. For 
the solutions shown, we find on the lower branch the global masses Mk=i ~ 2.357 and 
Mfc=3 ^ 2.371, and on the upper branch we find Mk=i ~ 10.25. Thus the fields outside 
the horizon contribute little to the mass. 

In Fig. 5b the local angular momentum J{x) is shown. Again, the global angular 
momentum, read off at infinity, is only slightly larger than the angular momentum 
at the horizon. For the solutions shown, the global angular momentum on the lower 
branch is Jk=i ~ 1.852 and Jk=3 ~ 1.872, and on the upper branch it is Jk=i ~ 102.5. 

The local electric charge Q{x) is shown in Fig. 5c. For solutions on both branches 
Q{x) decreases significantly with increasing x, giving rise to a large difference between 
the horizon electric charge and the global electric charge. In general Q{x) is smaller on 
the lower branch than on the upper branch, where it assumes values between ^ 0.9 at 
the horizon and 0.09 at infinity for the solutions shown. The decrease with x indicates 
the presence of a charge density outside the horizon, cancelling part of the horizon 
charge. 

We show the local magnetic charge P{x) in Fig. 5d. In contrast to the local electric 
charge, the local magnetic charge is not monotonic but possesses a maximum of magni- 
tude one, for the solutions on both branches. At infinity P{x) vanishes, in accordance 
with the boundary conditions, corresponding to a vanishing global magnetic charge. 
With increasing node number, the local magnetic charge remains close to the value one 
in an increasing region, as in the static case 0. This again indicates, that in the limit 
k —>■ oo, the rotating non-Abelian black hole solutions tend to a limiting solution with 
magnetic charge \P\ = 1, like their static counterparts 



5.5 Horizon properties 

We now turn to the horizon properties of the rotating non-Abelian black holes. In 
Figs. 6a-d we demonstrate the dependence of the horizon mass, the horizon angular 
momentum, and the horizon electric and magnetic charges on cjh for fixed isotropic 
horizon radius xh, for black hole solutions with one node. The corresponding horizon 
size, the deformation of the horizon and the surface gravity are shown in Figs. 7a-c. 

The horizon mass Ma follows closely the global mass M, as seen in Fig. 6a, where 
the horizon mass is shown for black hole solutions with horizon radius Xn = 0.1 and 
xh = 1. Similarly, the horizon angular momentum Ja closely follows the global an- 
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gular momentum J. The horizon angular momentum per unit mass, a a = Ja/M^, is 
shown in Fig. 6b. Both horizon mass and horizon angular momentum are close to the 
corresponding Kerr- Newman values for solutions with Q = and P = — 1. 

The horizon electric charge Q^, shown in Fig. 6c, increases monotonically along 
both branches, like the global non-Abelian electric charge Q. It is, however, about 
an order of magnitude larger than the global non-Abelian electric charge Q. The 
horizon magnetic charge Pa is shown in Fig. 6d. Whereas the black holes carry no 
global non-Abelian magnetic charge, their horizon magnetic charge is on the order of 
one. Starting from finite values on the lower branch, which correspond to the horizon 
magnetic charges of the static solutions. Pa decreases monotonically along the lower 
branch. It then reaches a minimum along the upper branch and increases again. The 
location of the minimum strongly depends on xh- 

The horizon size as quantified by the area parameter xa is shown in Fig. 7a for black 
holes with isotropic horizon radii xh = 0.1 and xh = 1 as a function of cuh- The horizon 
area grows monotonically along both branches, and diverges along the upper branch. 
Comparison with the corresponding Abelian black hole solutions again shows, that the 
non-Abelian horizon size follows closely the Kerr-Newman horizon size for Q = and 
P = — 1, whereas only for large xh it is also close to the Kerr horizon size. 

The deformation of the horizon is revealed by measuring the circumference of the 
horizon along the equator, Le, and the circumference of the horizon along the poles, 
Lp, Eqs. (|8lD. The deformation of the horizon, quantified by Le/Lp, Eq. (|8T[), is shown 
in Fig. 7b. The ratio Le/Lp grows monotonically along both branches. As tun tends to 
zero on the upper branch, the ratio tends to the value L^/ Lp 1.645, corresponding 
to the value of the limiting extremal Kerr solution, for M oo. On the lower branch 
the ratio L^/ Lp assumes the value one in the limit c<jh — > 0, corresponding to the value 
of a static spherically symmetric hairy black hole. Also shown are the corresponding 
embedded Abelian solutions, approaching the same limiting values of the deformation. 

The surface gravity Ksg is exhibited in Fig. 7c, along with the surface gravity of 
the Kerr and the Kerr-Newman solutions for Q = and P = — 1. The surface gravity 
of the non-Abelian black holes decreases monotonically along both branches, starting 
from the value of the corresponding static non-Abelian black hole solution on the lower 
branch in the limit cjh — >■ 0. On the upper branch the surface gravity tends to zero 
in the limit cuh — > 0, the value assumed by extremal black hole solutions. The surface 
gravity of Kerr-Newman and non-Abelian solutions agrees well for black hole solutions 
with large horizon radii, for black hole solutions with small horizon radii, however, a 
difference is observed on the lower branch. Here, in the limit cuh — >■ 0, the rotating 
non-Abelian black hole solutions approach the corresponding static non-Abelian black 
hole, and the Kerr-Newman solutions approach the corresponding Reissner-Nordstr0m 
black hole with Q = and P = — 1. For these static solutions the difference in surface 
gravity in known to decrease with increasing horizon size [Q. 
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5.6 Fixed cvn 



Having considered the rotating hairy black hole solutions for fixed horizon radius xu 
as a function of un, let us now keep cuh fixed and vary the horizon radius. 

In Fig. 8 we show the mass M of the non-Abelian black hole solutions as a function 
of the isotropic horizon radius xu for tun = 0.01, u^l = 0.02, and cuh = 0.05. For a given 
value of there is a minimal value of the horizon radius xr. In particular, the limit 
xh — is only reached for cur — >■ 0. Thus we do not obtain globally regular rotating 
solutions in the limit xh — ^ |]12| . 

For comparison. Fig. 8 also presents the mass M of the corresponding Kerr solutions 
and Kerr-Newman solutions with Q = and P = —1. For a fixed value of un, the 
mass of the Kerr solutions forms two straight lines, extending from the origin. The 
non-Abelian solutions tend toward these lines for large values of the horizon radius. 
The mass of the Kerr-Newman solutions again is close to the mass of the non-Abelian 
solutions. In particular, we observe, that the minimal values of the horizon radius only 
differ slightly for the non-Abelian and Kerr-Newman solutions. 



5.7 Fixed ujn/xn 

Let us finally consider variation of the parameters xh and cuh, while keeping their ratio 
cuh/xh fixed. 

The global charges of the black hole solutions with one node are shown in Fig. 9 for 
fixed ratio cjh/xh = 0.04. The mass M and the angular momentum per unit mass a 
increase monotonically along both branches, reaching finite limiting values on the upper 
branch for cuh — > 0. Except in a small region close to cuh = on the lower branch the 
non-Abelian electric charge Q also increases monotonically. The non-Abelian electric 
charge also approaches a finite limiting value on the upper branch, apparently close to 
the limiting value Qum ~ 0.124, observed previously. 

To gain a better understanding of the limiting behaviour on the upper branch, we 
consider also the horizon properties for these black hole solutions. In Fig. 10 we show 
their horizon size xa, their deformation L^/Lp and their surface gravity Kgg. On the 
upper branch the horizon size, the deformation and the surface gravity of the non- 
Abelian black hole solutions are very close to those of the Kerr-Newman solution with 
Q = and P = —1. In particular, the horizon size and deformation remain finite on 
the upper branch in the limit u^h ~^ 0, and the surface gravity tends to zero. This 
indicates, that an extremal black hole is approached in this limit. Since the limiting 
black hole solution retains its non-Abelian character in the gauge field functions, it 
should correspond to a rotating hairy extremal black hole. 

On the lower branch in the limit 0, the corresponding Bartnik-McKinnon 

solution is approached by the non-Abelian solutions, whereas the Kerr-Newman solu- 
tions approach a Reissner-Nordstr0m solution. This limiting behaviour is suggested by 
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a detailed inspection of the metric and gauge field functions. In Fig. 10 this limiting 
behaviour is reflected in the fact, that the horizon size of the non-Abelian solutions 
tends to zero, whereas the horizon size of the Kerr-Newman solutions tends to a fi- 
nite value. Furthermore, the surface gravity of the non-Abelian solutions apparently 
diverges in the limit, in agreement with the static non-Abelian results 0, whereas the 
surface gravity of the Kerr-Newman solutions tends to zero, the value of an extremal 
Reissner-Nordstr0m solution. 

6 Conclusions 

We have given a detailed account of a new class of black hole solutions in SU(2) EYM 
theory, which represent the first examples of non-perturbative stationary non-Abelian 
black hole solutions ||^. These black hole solutions carry mass, angular momentum 
and a non-Abelian electric charge. Although they do not carry a non-Abelian magnetic 
charge, they still possess non-trivial magnetic gauge fields outside their regular event 
horizon. They therefore represent rotating hairy black hole solutions. 

The global charges of the rotating hairy black hole solutions are not independent. 
For a given mass and angular momentum, as well as node number of the solution, 
a unique electric charge is obtained. Whereas mass and angular momentum are un- 
bounded, we observe, that the electric charge remains very small. 

The event horizon of the static axially symmetric black hole solutions resides at 
a surface of constant isotropic radial coordinate, x = x^. The boundary conditions 
at the horizon ensure regularity of the horizon. The horizon mass and horizon angu- 
lar momentum are only slightly smaller than the global mass and the global angular 
momentum of the black hole solutions. The horizon electric charge is, however, signifi- 
cantly larger than the global electric charge, and the solutions possess horizon magnetic 
charge of order one, whereas their global magnetic charge vanishes. 

The rotating hairy black hole solutions emerge from the static hairy black hole 
solutions in the limit of vanishing angular momentum. Since the static spherically 
symmetric black hole solutions form a sequence labelled by the node number k of 
the gauge field function, we obtain the corresponding sequence of rotating black hole 
solutions, by starting from the static black hole solutions and imposing a small angular 
velocity of the horizon via the boundary conditions. 

The rotating non-Abelian black hole solutions are rather close to the Kerr-Newman 
solutions with Q = and P = —1. In particular, the metric and gauge field functions 
of the non-Abelian solutions are very close to those of the Abelian solutions, except for 
those gauge field function, which do not vanish in the static limit. These retain the 
non-Abelian character of the solutions. However, with increasing node number, also 
these functions tend to their Abelian counterparts. 
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The asymptotic expansion performed for the metric and gauge field functions, con- 
tains non-integer powers for the magnetic gauge field functions. In particular, the 
non-integer exponents depend on the non-Abelian electric charge. The expansion then 
imposes constraints on the possible values of the electric charge. Since in the static 
limit the non-Abelian electric charge vanishes, the well known power law decay of the 
static gauge field functions is recovered. 

The expansions of the metric and gauge field functions at the horizon show, that the 



rotating hairy black hole solutions satisfy the zeroth law of black hole mechanics |19 



Relations obtained recently within the isolated horizon framework |T^, ^ concerning 
various horizon properties, such as the horizon mass and the horizon charges, will be 
considered elsewhere. 

The hairy black hole solutions constructed here non-perturbatively, were first con- 
sidered perturbatively [1^. To compare with these perturbative calculations, where 



linear rotational excitations of the static EYM black holes were studied, we consider 
the slowly rotating non-Abelian solutions in the limit cjh — > 0. In the perturbative 
calculations Q oc J [|l^, and the ratio Q/J depends only on the horizon radius. The 



non-perturbative calculations show good agreement with the non-perturbative results 
for the slowly rotating solutions with large values of the horizon radius. 

Besides these non-Abelian stationary black hole solutions with finite angular mo- 
mentum J and finite electric charge Q, perturbative studies |]l5l have predicted two 
more types of stationary non-Abelian black hole solutions. These correspond to rotat- 
ing black hole solutions which are uncharged (J > 0, Q = 0), and non-static charged 
black hole solutions, which have vanishing angular momentum (J = 0, Q 7^ 0). Both 
types satisfy a different set of boundary conditions at infinity. 

This different set of boundary conditions at infinity should also be observed by 
rotating regular non-Abelian solutions [|15|. The numerical construction of such non- 



perturbative regular solutions has been attempted recently [28|, and arguments have 



been put forward, that such solutions should not exist. Our attempts to obtain numer- 
ically the non-perturbative counterparts of the predicted further types of black hole 
solutions have met with the same difficulties for the same reasons 

Besides the rotating black hole solutions considered here, there might be rapidly 
rotating branches of non-Abelian black holes solutions in EYM theory, not connected 
to the static solutions. There should also be rotating hairy black hole solutions in 
other non-Abelian theories such as Einstein- Yang-Mills-Higgs theory. For instance, we 
expect rotating black hole solutions with magnetic charge and with magnetic dipole 
hair 0, ^, Furthermore, the recent conjecture, that "any dyon solution with nonzero 
angular momentum necessarily contains an event horizon" p8[ awaits investigation. 
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7 Appendix A: Ricci circularity conditions 



Let us demonstrate that the ansatz for the metric @ satisfies the Ricci circularity 
conditions, Eq. ([Tl|) . 

The Kilhng vectors C, = dt and t] = have components 

e = (1, 0, 0, 0) , = C9^.u = {gtt, 0, 0, gt^) , 

r = (0, 0, 0, 1) , v, = v'g,. = (g^t, o, o, g^^) . (ii7) 

Consider the Ricci circularity condition i^R^ia^pV-y] = 0; Eq. (|ll]), where the square 
bracket denotes antisymmetrization. Hence, a ^ (i ^ ^. Thus i^R^iaip'U'^] 7^ is only 
possible if either /3 = 0, 7 = 3 or /5 = 3, 7 = 0. In both cases a = 1 or a = 2. 
For a = 1, 

^^R^l[liQV2.] = Ro[i^oV3] 

= ^{Roi^oVs + Roo^sVi + Ros^iVo - Ros^oVi - Roo^m - i?oi6^o) 


= 2^01(^0^3-6^0) 

= -^Roiigugipv ~ gl,t) ^ 

and similarly, for a = 2, 

^''Ri,i2^oV3] = ^Ro2{gttgw - alt) , 

since all other components vanish due to6 = 0,6 = and ?7i = 0, 772 = 0. 
On the other hand, the ansatz for the metric Eq. (|]) yields 

Rqi = Rtr = , R02 = Rte = . 

Consequently, ^^R^ia^pi]^] = 0. 

Note that Rtr = 0, Rto = implies Gtr = Rtr — ^f^gtrR = and Gte = Rte — 
l/2gteR = 0, respectively. Therefore Ttr = and Tte = is a necessary condition for 
the solutions of the Einstein equations. In a similar fashion it can be shown that the 
Ricci circularity condition i]^ R^\^a^prj-^^ = 0, Eq. ([Tl|) , is satisfied by the metric (§) and 
implies T^r = and T^e = 0. However, for the Ansatz ([T5|) for the gauge field these 
conditions are fuUfiUed identically. Hence the solutions satisfy the Ricci circularity 
conditions. 
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8 Appendix B: Asymptotic expansion 



8.1 General asymptotic expansion 



The asymptotic expansion of the fields may be obtained from the field equations and 
the corresponding boundary conditions. However, the process is rather involved. In- 
deed, the most natural assumption for the asymptotic r-dependence of the functions, 
i.e. polynomial, seems to be in contradiction with the presence of a non-vanishing elec- 
tric charge. In view of that and inspired by perturbative results fl^, we allowed for the 
presence of logarithms in the expansions. We then observed that, when keeping only 
a finite number of logarithmic terms, the charge was forced to vanish. By permitting 
an infinite number of such terms, however, it became possible to obtain a consistent 
expansion in the presence of an electric charge. This feature then suggested to include 
non-integer powers of r, with the exponents depending on the electric charge. 

Numerically, we realized that there should be a splitting of some terms in the 1/r 
Taylor series expansion of the static spherically symmetric case. However, the exponents 
of the new terms turned out to be very close to integer numbers, even though the 
behavior of the functions could not be described just by means of a 1/r Taylor series. 
The procedure of how to compute this expansion was then clear. 

First of all, we introduced a formal parameter e in order to characterize terms of 
the order of 1/r, without assuming a Taylor series in 1/r for the functions. Then we 
expanded all the functions in this formal parameter, the coefficients of such expansions 
depending on r and 9. Finally, we introduced those e series into the system of field 
equations, taking into account the explicit dependence on r in the equations by includ- 
ing e appropriately. The last step was collecting coefficients in e for each equation of the 
system, and solving the equations so formed, order by order in this formal parameter, 
keeping in mind the boundary conditions and the fact that the coefficients of the e 
series had to behave consistently with the corresponding power of e. 

Proceeding in this way, the expression for this asymptotic expansion is found to be: 
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with the notation of section 3.1. 
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8.2 Asymptotic expansion: relation with static case 

Here we show how the previous asymptotic expansions reduce to the ones correspond- 
ing to the static case in the hmit of vanishing charge. Moreover, as these are non- 
perturbative expansions, they must include the perturbative expansions reported in 
Jl^. This is indeed the case. In order to perform perturbative expansions in terms 



of the electric charge, we must recall that, for a fixed value of the horizon radius, all 
the parameters in these expansions are functions of Q. Due to the presence of Q in 
some denominators, the dependence on Q of these coefficients has to be singular, in 
such a way that the resulting series in Q turns out to be regular. The behavior of the 
constants as functions of Q reads: 

Mo + Q2Ko(g) , 

-]^Ml + Q^K,{Q) , 

-h + Q^K2{Q) , 
h + Q^K,{Q) , 

4Q2 



M 

C2 
C4 



C5 



1 52 

1 62 
QK,{Q) , 



h,2 + ih{K,{Q)-K,m-\Km 



+ QMQ) 



16 



4 



20 



(119) 

where Mq is the dimensionless mass of the static spherically symmetric solution, b is 
the parameter of the gauge field for such a limiting solution, and i^i, . . . , are regular 
functions of Q. Introducing these relations in ( |118|) and expanding the result in Q, we 
recover the perturbative expansions given in UTT 
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where primes denote the notation of Volkov and Straumann, which is related to our 
notation by 
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8.3 Asymptotic expansion: embedded Kerr-Newman solutions 

For comparison, we present the asymptotic expansion also for the embedded Kerr- 
Newman solutions for our choice of coordinates and gauge. The expansion reads 
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2aM sin 6* ^ / 1 , , , 

B2 = + , (129) 



where a is the angular momentum per unit mass, M is the mass, Q is the electric 
charge, (Q^ = Q"'Q°'), and P is the magnetic charge, (P^ = P"-P"-). 



9 Appendix C: Expansion at the horizon 

We here first motivate our choice of boundary conditions at the horizon. Then we give 
the full expansion of the metric and gauge field functions at the horizon and relate this 
general expansion with the static case. 



9.1 Boundary conditions at the horizon 



Let us begin by noting that the ansatz of the gauge field has the property [|T^, |2^ 

d^A^ = D^u , (130) 
with u = Tz/2. The components F^^ can be expressed as 



with 



W = A^-u 



(131) 



(132) 



W transforms as a scalar doublet under gauge transformations U = exp(irr<^/2), 
Eq. (pOl) . Using the definition of W, Eq. (|132|) , we find for the component At of 
the gauge field 

Ij) T I,} ^ 1,1 

(133) 



A = <ij + -^ + -w = i' + -w 

r 2 r r 



with 



^ = ^ + ^^-^ 
r 2 



(134) 



Thus also transforms as a scalar doublet under gauge transformations. 

To discuss the behaviour of the solutions at the horizon it is convenient to rewrite 



^ and A^ as 
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This yields for the scalar doublets, W and Eqs. (|132|) , (|134|) , 
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We assume (see below) that near the horizon the metric functions can be expanded 



as 



:i38i 



where 6 = {x — xii)/xu and is a constant. For the gauge field functions we assume 
an expansion in the form 
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(139) 



where -ffi|a;=XH = fi^^s the gauge freedom. 
Let us write the field equations as 



-9 



-gF^n = , 



fJ,U 
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The expansion of E^ at the horizon yields 

Bu = Hio{uJu -uji)sm9 , 
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With this result the expansion of -E*^ leads to the conditions 
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In terms of W and ^ these conditions are equivalent to 
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We now assume that the electro-static potential \1/ is constant at the horizon, 
Eq. (|55D , i. e. Biq = const and B20 = 0. 

To discuss the boundary conditions ( |143| ) , let us first assume '^1^=^^ ^ and 
l^^l^j^ajj^ = \^\x=XYii for some function \{ff). In this case the expansion yields for the 
gauge potential 
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The functions Hi, 1 — H2, H3 and B2 vanish at least up to order 0{6^), indicating that 
we find only embedded Abelian solutions for ^[^.^j-jj ^ 0. 

Assuming "i/l^^^^ = 0, however, does not imply restrictions on H20, H^o and H^q. In 
this case the expansion yields H21 — H31 — H41 — and coi — uju- The last condition 
implies Bu — B21 — 0. In this case non-Abelian solutions are possible. 



9.2 General expansion at the horizon 

Here we present the expansion of the functions of the stationary axially symmetric 
black hole solutions at the horizon in powers of 6. These expansions can be obtained 
from the regularity conditions imposed on the Einstein equations and the matter field 
equations: 
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B12 and i?22 are functions of 6. Relations ([76|) and ( [77D also hold. 
9.3 Expansion: relation with static case 

The general expansion at the horizon includes the expansion for the static case. The 
static limit corresponds to setting = uj2 = -B12 = -B22 = -f^ii = -f^ai = 0, m2 = h, 
and H4Q = H20] in addition, /2, I2, and H20 become constant. The expansion reads 
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recovering the known expressions for the behavior at the horizon of the static spherically 
symmetric black hole solutions. 



9.4 Expansion: embedded Kerr-Newman solutions 

For comparison, we present the expansion at the horizon also for the embedded Kerr- 
Newman solutions for our choice of coordinates and gauge. The expansion reads 
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^^h[(^ + 2.h)^ + «^-«^^]^^2(i_^)^0(0, (147) 



[2M(M + 2xh) - + P^)Y 



4[(M + 2.H)^ + a^cos^.p^^,^^_3^^^^^^,^^ (148) 



[2M(M + 2a;H) - (Q^ + p2)]2 
45^(1 -35) + 0(5^) , (149) 

(1 + 5) 



2M(M + 2a;H) - (Q^ + p2) 
2aa;?r 



-{(M + 2xh)[2(M + 2xh)(M + xn) 



[2M(M + 2a;H) - (Q^ + p2)]3 

-(g' + P2)] + 2a2a;H cos' + 0{5^) , (150) 

, (151) 

, (152) 

- (M + 2.H)' + a^cos^^ { ' + + ^[2^^^ + 

-(g2 + P2)] + a2 ^Qg2 ^ ^g^^ ^ 2xh) sin ^} 

-[2a'gxH cos' - 4aPxH(M + 2xh) cos^ 



2[(M + 2xH)' + a2cos2 
-2Qxh{M + 2xh)']5' + 0(5^) , (153) 
, (154) 

g(M + 2xH) - a cos 6' 
2M(M + 2xn) - (Q' + P') 

- 2|2M(M + 2.H)-(Q- + P')p {^'^-" [ - '^<^ ^ '^«' - ^ ^'^'^ 

+ (M + 2a;H)'[(M + 2xnf + 44]] - Aaxii[{M + 2a;H)[4MPxH 

-(g' + P2)] - 2PxH(g' + P') + 2(M + 2a;H)'(M + xn)] cos ^ 

+8a'g4(^ + 2xh) cos' ^ - 8a^4 cos=^ ^}<^^ + 0(5^) , (155) 

(I sill ^? 

2M(M + 2xn) - {Q^ + P') 

9<TrTT sin or, 

r{(M + 4a;H)[2M(M + 2a;H) - {Q^ + P')] 



[2M(M + 2xn) - {Q^ + P')]^ 
-2a'xH sin' ^}5' + 0{6^) , (156) 



= ^V^'-(a2 + g2 + p2) , 
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(157) 



S = ^-1 



(158) 



42 



References 



[1] W. Israel, Commun. Math. Phys. 8 (1968) 245; 
D. C. Robinson, Phys. Rev. Lett. 34 (1975) 905; 
P. Mazur, J. Phys. A15 (1982) 3173. 

[2] M. Heusler, Black Hole Uniqueness Theorems, (Cambrigde University Press, 
1996). 

[3] see e.g. M. S. Volkov and D. V. Gal'tsov, Phys. Rept. 319 (1999) 1; 

D. V. Gal'tsov, [hep-th/0 112038] , Proceedings of the 16th International Conference 
on General Relativity and Gravitation, July 2001, Durban, South Africa 

[4] M. S. Volkov and D. V. Galt'sov, Sov. J. Nucl. Phys. 51 (1990) 747; 
P. Bizon, Phys. Rev. Lett. 64 (1990) 2844; 

H. P. Kiinzle and A. K. M. Masoud-ul-Alam, J. Math. Phys. 31 (1990) 928. 

[5] B. Kleihaus and J. Kunz, Phys. Rev. Lett. 79 (1997) 1595; 
B. Kleihaus and J. Kunz, Phys. Rev. D57 (1998) 6138. 

[6] B. Kleihaus and J. Kunz, Phys. Rev. Lett. 85 (2000) 2430; 
B. Kleihaus and J. Kunz, Phys. Lett. B494 (2000) 130. 

[7] B. Hartmann, B. Kleihaus, and J. Kunz, Phys. Rev. Lett. 86, 1423 (2001); 
B. Hartmann, B. Kleihaus, and J. Kunz, Phys. Rev. D65 (2002) 0024027. 

[8] D. Sudarsky and R. M. Wald, Phys. Rev. D46 (1992) 1453; 
D. Sudarsky and R. M. Wald, Phys. Rev. D47 (1993) 5209. 

[9] see e. g. D. Kramer, H. Stephani, E. Herlt, and M. MacCallum, Exact Solutions 
of Einstein's Field Equations, Ch. 17 (Cambridge University Press, Cambridge, 
1980) 

[10] A. Papapetrou, Ann. Inst. H. Poincare A4 1966 83; 
B. Carter, J. Math. Phys. 10 (1969) 70. 

[11] M. Heusler and N. Straumann, Class. Quant. Grav. 10, (1993) 21. 

[12] M. S. Volkov and N. Straumann, Phys. Rev. Lett. 79 (1997) 1428. 

[13] D. V. Gal'tsov, Einstein- Yang-Mills solitons: towards new degrees of freedom, 
|gr-qc/980gim . 

[14] D. V. Gal'tsov and A. A. Ershov, Phys. Lett. 138 A (1989) 160. 



43 



[15] O. Brodbeck, M. Heusler, N. Straumann and M. S. Volkov, Phys. Rev. Lett. 79 
(1997) 4310; 

O. Brodbeck and M. Heusler, Phys. Rev. D56 (1997) 6278. 
[16] B. Kleihaus and J. Kunz, Phys. Rev. Lett. 86 (2001) 3704. 

[17] A. Corichi, and D. Sudarsky, Phys. Rev. D61 (2000) 101501; 

A. Corichi, U. Nucamendi, and D. Sudarsky, Phys. Rev. D62 (2000) 044046. 

[18] B. Kleihaus and J. Kunz, Phys. Rev. Lett. 78 (1997) 2527; 

B. Kleihaus and J. Kunz, Phys. Rev. D57 (1998) 834. 

[19] R. M. Wald, General Relativity (University of Chicago Press, Chicago, 1984) 

[20] J. N. Islam, Rotating fields in general relativity (Cambridge University Press, 
Cambridge, 1985) 

[21] A. Ashtekar, Ch. Beetle and J. Lewandowski, Phys. Rev. D64 (2001) 044016. 
[22] M. J. Perry, Phys. Lett. B71 (1977) 234. 

[23] When both electric and magnetic non-Abelian charge is present, the vectors must 



be parallel to each other, i.e. Q"" = AP" with A constant 



[24] Since the gauge fields wind n times around, while the azimuthal angle if covers 
the full trigonometric circle once, we refer to the integer n as the winding number 
of the solutions. 

[25] B. Kleihaus, J. Kunz, A. Sood and M. Wirschins, Phys. Rev. D58 (1998) 084006. 

[26] The limiting solution corresponds to an embedded Reissner-Nordstr0m black hole 
with magnetic charge \P\ = n for horizon radius xh > 1, where x represents a 
Schwarzschild-like radial coordinate. [0, |^. For xh < 1, the limiting solution also 
corresponds to an embedded Reissner-Nordstr0m black hole with magnetic charge 
\P\ = n in the region 1 < x, but contains an oscillating non-Abelian part in the 
region xh < x < 1 [0, 

[27] W. Schdnauer and R. WeiB, J. Comput. Appl. Math. 27, 279 (1989) 279; M. 
Schauder, R. WeiB and W. Schonauer, The CADSOL Program Package, Univer- 
sitat Karlsruhe, Interner Bericht Nr. 46/92 (1992). 

[28] J. J. van der Bij and E. Radu, On rotating regular nonabelian solutions. Int. J. 
Mod. Phys. A17 (2002) 1477. 

[29] B. Kleihaus, J. Kunz, F. Navarro, in preparation. 



44 



Fig. la 




Figure la: The dimensionless mass M is shown as a function of oj-a for node number 
k — 1 and xu — 1.0, 0.5 and 0.1. For the same values of parameters the dimensionless 
mass of the Kerr solution (thin solid) and the Kerr-Newman solution (dotted) for Q — 
and \P\ = 1 are also shown. 
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Fig. lb 




Figure lb: Same as Fig. la for the ratio a — J/M. 

Fig. Ic 




Figure Ic: The electric charge Q is shown as a function of uj-r for node number k — 1 
and xh = 1.0, 0.5 and 0.1. 
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Fig.2 




Fig. 3 
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Figure 3: The same as Fig. 2 on the upper branch. 
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Fig. 4a 
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Figure 4a: The function H\ is shown for = 1-0, = 0.04 and A; = 1 on the lower 
and upper branch, and for /c = 3 on the lower branch. 

Fig. 4b 
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Figure 4b: Same as Fig. 4a for the function i?2- 
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Fig. 4c 
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Figure 4c: Same as Fig. 4a for the function H3. The curves for A; = 1 and A; = 3 on 
the lower branch coincide. 

Fig. 4d 
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Figure 4d: Same as Fig. 5a for the function H4. 
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Fig. 4e 
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Figure 4e: Same as Fig. 4a for the function Bi. 
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Figure 4f: Same as Fig. 4a for the function B2. 
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Fig. 4g 
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Figure 4g: Same as Fig. 4a for the function /. The curves for A; = 1 and A; = 3 on the 
lower branch coincide. 

Fig. 4h 
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Figure 4h: Same as Fig. 4a for the function g = m/l. The curves for k = 1 and /c = 3 
on the lower branch coincide. 
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Fig. 4i 
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Figure 4i: Same as Fig. 4a for the function /. The curves for k — 1 and A; = 3 on the 
lower branch and A; = 1 on the upper branch coincide. 

Fig. 4j 
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Figure 4j: Same as Fig. 4a for the function a;. The curves for k = 1 and /c = 3 on the 
lower branch coincide. 
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Fig. 5a 
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Figure 5a: The local mass M{x) is shown for xu = 1.0, a;H = 0.04 and A; = 1 on the 
lower and upper branch, and for /c = 3 on the lower branch. 

Fig. 5b 
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Figure 5b: The same as Fig. 5a for the local angular monetum J{x). 
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Fig. 5c 
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Figure 5c: The same as Fig. 5a for the local electric charge Q{x). 



Fig. 5d 




Figure 5d: The same as Fig. 5a for the local magnetic charge P{x). 
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Fig. 6a 
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Figure 6a: The horizon mass Ma is shown as a function oi cuu ior k — 1, xn — i and 
xn — 0.1 on the lower branch (solid) and on the upper branch (dashed). For the same 
values of parameters the corresponding functions of the Kerr solution (thin solid) and 
the Kerr-Newman solution (dotted) for Q = and |P| = 1 are also shown. 
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Fig. 6b 
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Figure 6b: The same as Fig. 6a for the angular momentum per mass at the horizon 
a A = J a/ Ma- 

Fig. 6c 
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Figure 6c: The same as Fig. 6a for the horizon electric charge Qa- 
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Fig. 6d 
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Figure 6d: The same as Fig. 6a for the horizon magnetic charge Pa- 
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Fig. 7a 




Figure 7a: The area parameter is shown as a function of a;H for /c = 1, xh = 1 and 
xh = 0.1 on the lower branch (sohd) and on the upper branch (dashed). For the same 
values of parameters the corresponding functions of the Kerr solution (thin solid) and 
the Kerr-Newman solution (dotted) for Q = and |P| = 1 are also shown. 
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Fig. 7b 




Figure 7b: The same as Fig. 7a for the horizon electric charge Qa- 



Fig. 7c 




Figure 7c: The same as Fig. 7a for the surface gravity n. 
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Fig. 8 
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Figure 8: The dimensionless mass M is shown as a function of xu ior k — 1 and 
LOfi — 0.01, 0.02 and 0.05 on the lower branch (sohd) and on the upper branch (dashed). 
For the same values of parameters the corresponding functions of the Kerr solution (thin 
solid) the Kerr- Newman solution (dotted) for Q = and |P| = 1 are also shown. 
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Fig. 9 




Figure 9: The dimensionless mass M , ratio a = J/M and electric charge Q are shown 
as function of for k = 1 and fixed uJn/xn = 0.04 on the lower branch (solid) and 
on the upper branch (dashed). For the same values of parameters the corresponding 
functions of the Kerr solution (thin sohd) the Kerr- Newman solution (dotted) for Q = 
and IPI = 1 are also shown. 
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Fig. 10 




Figure 10: The same as Fig. 9 for the horizon parameter x^, the ratio Le/Lp and the 
surface curvature k. 
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